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ABSTRACT 


The  magnetohydrodynamic  equations  for  temporal 
transients  have  been  formulated  and  solved  for  a  liquid 
metal  flowing  in  a  rectangular  channel.  The  rectangular 
channel  has  a  perfectly  conducting  moving  top  wall  and  a 
perfectly  conducting  stationary  bottom  wall  in  the  presence 
of  an  applied  external  magnetic  field  aligned  perpendicular 
to  the  conducting  walls.  The  side  walls  of  the  channel  are 
stationary  insulators.  Calculations  show  that  the  temporal 
transients  of  the  fluid  velocity  and  induced  magnetic 
intensity  comprise  two  exponentially  decaying  parts.  One  fast 
transient  is  believed  to  be  associated  with  the  propagation  of 
Alven  waves  and  the  other  a  slow  transient  being  the  result 
of  viscous  and  electrical  diffusion.  Curves  of  the  transients 
are  presented  at  several  stations  in  the  channel. 

ADMINISTRATIVE  INFORMATION 

This  work  was  a  cooperative  effort  between  the  David  Taylor  Research 
Center  and  the  Scientific  Division  of  Frontier  Timber  Co.,  Bradford,  Pa.  It  was 
performed  under  Program  Element  63224C,  Work  Unit  1-2712-501,  Project 
Title:  SDI  Pulse  Power  Key  Technology,  Responsible  Individual  P.  Filios.  The 
work  was  also  partially  supported  by  the  DTRC  Independent  Research  Program, 
Director  of  Naval  Research,  OCNRIO,  and  administered  by  the  Research 
Director,  DTRC0113  under  Program  Element  61152N,  Project  Number 
ZR00001,  Task  Area  ZR0230201,  Work  Unit  1-2712-125,  Project  Title: 
Orientation  Effects  in  Liquid-Metal  Collectors. 

INTRODUCTION 

The  use  of  liquid  metals  for  current  collectors  in  homopolar  motors  and 
generators  has  led  to  the  design  of  machines  of  superior  performance.  The 
steady  state  power  losses  have  been  studied  in  a  model  of  a  current  collector 
comprising  a  liquid-metal  rectangular  channel  with  an  applied  external  magnetic 
field  and  boundary  conditions  containing  a  combination  of  moving  and  fixed 
insulating  and  perfectly  conducting  walls.1  The  top  moving  wall  and  stationary 
bottom  wall  are  perfect  conductors.  The  side  walls  are  insulators.  The  applied 
external  magnetic  field  is  perpendicular  to  the  conducting  walls. 

In  some  applications  of  homopolar  generators  it  becomes  necessary  not 
only  to  start  and  stop  the  machines  but  also  to  operate  them  under  oscillating 
conditions.  This  could  be  the  case  in  an  application  where  a  homopolar 
generator  behaves  as  an  extremely  high  energy  capacitor.  Therefore,  one  is 
interested  in  examining  energy  losses  caused  by  these  non-steady  state 
phenomena. 

Studies  of  non-steady  state  Couette  flow  in  an  external  magnetic  field  are 
seldom  found  in  the  literature.  There  have  been  some  papers  which  augmented 
the  classical  Rayleigh  problem2  with  an  external  magnetic  field  in  an  electrically 
conducting  medium.3-8  The  presence  of  Alven  waves  and  slowly  propagating 
waves  was  observed  in  detailed  studies  of  the  Rayleigh  problem.  Theory 
predicted  the  possible  formation  of  switch-on  shocks  and  current  sheets  to  be 


present  in  the  Alven  wave  front  under  certain  conditions.  Currently,  it  would 
appear  that  a  complete  theoretical  analysis  and  explanation  of  the  various 
phenomena  comprising  the  Rayleigh  problem  is  lacking.  This  is  partly  caused  by 
a  waning  interest  in  the  field  of  magnetohydrodynamics  in  the  past  decades  and 
the  complexity  of  the  mathematical  solutions  in  these  magnetohydrodynamic 
problems. 

Our  objective  is  to  present  a  transient  two-dimensional  solution  and  some 
numerical  results  for  one  particular  geometry  of  finite  extent,  shown  in  Fig.  1. 
The  solution  enables  the  investigation  of  the  influence  of  the  sidewalls  on  the 
velocity  and  the  induced  magnetic  field  that  affects  the  joulean  and  viscous 
dissipation.  The  geometry  chosen  roughly  corresponds  to  an  axisymmetric  liquid 
metal  current  collector  having  channel  dimensions  that  are  small  compared  to 
the  radius  of  curvature  of  the  current  collector.  The  moving  and  fixed 
conducting  walls  correspond  to  the  rotor  and  the  stator,  respectively.  Following 
earlier  work,1  the  orientation  of  the  applied  magnetic  field  or  the  applied 
induction,  B0  (see  Fig.  1)  is  chosen  to  generate  the  maximum  eddy  current.  In 
addition,  current  parallel  to  the  applied  magnetic  field  produces  no 
electromagnetic  body  force  and  need  not  be  considered.  Analytical  expressions 
are  derived  for  the  induced  magnetic  field  and  the  fluid  velocity  in  terms  of 
their  spatial  and  temporal  variations.  The  parameters  of  variation  are  applied 
magnetic  induction,  fluid  viscosity,  density  and  electrical  conductivity,  and  the 
channel  aspect  ratio.  Numerical  results  showing  the  buildup  of  the  velocity  and 
induced  magnetic  field  transients  are  presented  for  several  cases. 


GENERAL  TRANSIENT  MAGNETOHYDRODYNAMIC  EQUATIONS 

Consider  the  case  of  laminar  flow  of  an  incompressible,  uniformly 
conducting  fluid  of  density  q,  viscosity  /if,  and  conductivity  o  flowing  in  a 
rectangular  duct  of  width  2zq  and  height  of  y0  with  a  uniform  external  magnetic 
induction,  B0  applied  in  the  y  direction  as  shown  in  Fig.  1.  The  flow  is  driven 
by  the  sudden  movement  of  the  top  conductor.  This  conductor,  initially  at  rest, 
is  assumed  to  attain  a  constant  velocity  at  the  onset  of  the  transient  calculation 
and  maintains  that  velocity.  It  is  assumed  that  the  pressure  gradient  in  the  fluid 
is  zero  to  simplify  the  solutions.  Such  would  not  be  the  case  if  there  were  any 
obstructions  in  the  channel  that  limited  the  net  flow  of  fluid.  It  is  assumed  that 
no  secondary  flow  is  generated  during  the  transient,  and  there  is  no  variation  in 
the  duct  cross  section  or  distortion  of  the  applied  induction  in  accordance  with 
the  assumption  of  a  fluid  velocity  and  an  induced  magnetic  field  intensity  in  the 
x  direction  only.  The  equation  for  the  fluid  velocity  in  the  x  direction  is  given 
by 


6  3U*  _  32U*  32U*  3 H*  3 P* 

3t*  ^  3y*2  +  ^  dz*^  +  o  3 y*  3x* 


and  the  induced  magnetic  intensity  in  the  x  direction  is  given  by  the  magnetic 
diffusion  equation  as 


3H*  32H*  32H*  _  3U* 

OMo  3t*  -  3y*2  +  3z*2  +  °B°  3y* 

2 


(2) 


where  all  dimensional  variables  are  denoted  with  a  superscript  asterisk,  contrary 
to  earlier  work.1-8  Equation  (1)  is  the  two-dimensional  Navier-Stokes  equation 
including  an  electromagnetic  body  force  and  Eq.  (2)  is  the  magnetic  diffusion 
equation  obtained  by  combining  both  of  Maxwell’s  curl  equations  and  the 
electrical  constitutive  equation  for  a  moving  conducting  medium.  All  equations 
are  in  the  SI  or  RMKSA  system  of  units. 

RECTANGULAR  CHANNEL  WITH  MOVING  WALL 

The  problem  solved  in  this  paper  considers  a  rectangular  channel  (see 
Fig.  1)  filled  with  a  liquid  metal  in  a  uniform  magnetic  field,  B0.  The  top  wall  is 
assumed  to  be  perfectly  conducting  and  is  initially  at  rest.  At  time  t*  =  0  its 
velocity  jumps  in  a  unit  step  to  U0,  where  it  remains  constant.  The  bottom  wall 
is  stationary  and  is  also  an  ideal  conductor  while  the  side  walls  parallel  to  the 
applied  external  magnetic  induction  field  are  stationary  electrical  insulators.  To 
make  the  handling  of  the  equations  easier  we  define  a  number  of  convenient 
and  physically  meaningful  dimensionless  quantities.  The  dimensionless  quantities 
are  not  denoted  with  asterisks  as  in  previous  work  and  are  defined  as 

U  =  U*/U0,  y  =  y*/y0,  z  =  z*/y0,  t  =  t*/(y 0V/if), 

H  =  and  k  =  z0/yo»  the  Hartmann  number 

M  =  y0B0\/  a/ix{  ,  G  =  [ y02/(MfU0)]3P*/dx * 

the  magnetic  Prandtl  number  is  given  by  Pm  =  o^o^f/e,  where  jig  is  the 
permeability  of  a  vacuum.  Even  though  the  Hartmann  number  was  undefined  in 
the  magnetohydrodynamic  Raleigh  problem  (because  y0  was  unbounded),  it  is 
the  most  important  parameter  in  this  problem.  It  is  a  measure  of  the  ratio  of 
ponderomotive  force  to  the  viscous  force  (in  this  interpretation  it  is  assumed  the 
ponderomotive  force  is  of  the  order  of  U0,  certainly  true  for  the  electrical 
conductivities  encountered  in  liquid  metal  current  collectors).  The  magnetic 
Prandtl  number  is  small  in  this  problem  but  cannot  be  neglected  entirely  without 
causing  great  inaccuracy  in  the  solutions.  It  can  be  discarded  only  in  direct 
comparison  with  other  larger  quantities.  The  magnetic  Prandtl  number  is  a 
measure  of  the  ratio  of  vorticity  diffusion  to  magnetic  diffusion.  For  liquid  NaK 
the  magnetic  Prandtl  number  is  of  the  order  of  about  2x  10-6.  If  liquid  metal 
electrical  superconductors  were  available,  the  magnetic  Prandtl  number  would  be 
very  large  and  could  cause  interesting  effects.  To  make  the  solutions  much 
simpler  we  assume  G  =  0.  Then  Eqs.  (1)  and  (2)  become 

32U/3y2  +  32U/3z2  +  M3H/3y  =  3U/3t  (3) 

3 2H/ 3 y2  +  32H/3z2  +  M3U/3y  =  Pm  3H/3t  (4) 

Here  we  note  that  the  magnetic  Prandtl  number  provides  a  coupling  between  the 
induced  magnetic  intensity  and  the  fluid  velocity  that  does  not  exist  in  the 
steady  problem  treated  earlier.1  Equations  (3)  and  (4)  must  be  solved 
simultaneously  for  U(y,z,t)  and  H(y,z,t)  by  any  method  which  satisfies  the 


boundary  and  initial  conditions  while  satisfying  Eqs.  (3)  and  (4).  It  is  tempting 
to  consider  the  use  of  the  Laplace  transform  to  eliminate  time  as  a  variable  until 
the  difficulty  of  the  inversion  of  the  transform  is  considered.  Instead  of  using 
the  Laplace  transform  method  for  solution  of  the  partial  differential  equations, 
we  have  decided  to  seek  a  solution  in  terms  of  exponentially  diminishing  time 
functions.  To  do  that  we  let 

U(y,z,t)  =  Uss(y,z)  +  u(y,z,t)  (5) 

H(y,z,t)  =  Hss(y,z)  +  h(y,z,t)  (6) 

where  Uss  and  Hss  are  the  solutions  to  the  steady  flow  problem  and  u(y,z,t)  and 
h(y,z,t)  are  transients  that  eventually  fade.  When  Eqs.  (5)  and  (6)  are  substituted 
into  Eqs.  (3)  and  (4)  there  results  for  the  transients  u  and  h 

32u/3y2  +  32u/3z2  +  M3h/3y  =  3u/3t  (7) 

32h/3y2  +  32h/3z2  +  M3u/3y  =  Pm  3h/3t  (8) 

BOUNDARY  CONDITIONS 

The  no  slip  condition  applies  to  the  total  solution  for  the  velocity  at  all  of 
the  boundaries.  In  that  case 

U(y=  l,z,t)  =  1  implying  that  u(y  =  1  ,z,t)  =0 

U(y  =  0,z,t)  =  0  implying  that  u(y  =  0,z,t)  =  0 

U(y,z=  ±k,t)  =  0  implying  that  u(y,z=  ±k,t)  =  0 

Because  there  must  be  zero  electric  field  in  the  rest  frame  of  ,he  perfect 

conductors  for  finite  current  flow,  the  normal  derivative  of  the  magnetic 
intensity  is  zero  at  the  top  and  bottom  perfect  conductors.  Currents  generated  in 
the  top  conductor  and  the  fluid  move  in  a  loop  (in  the  y-z  plane)  bounded  by 
the  conductors  and  the  side  insulating  walls.  All  the  current  generated  by 
conductors  moving  in  the  applied  external  induction  field  returns  primarily  in 
the  bottom  perfect  conductor.  Thus  by  Ampere’s  law  the  induced  magnetic 
intensity  at  the  insulating  walls  is  zero.  In  that  case 

3H/3y|y=1  =  0  implying  that  3h/3y|y=1  =  0 

3H/3y|y  =  0  =  0  implying  that  3h/3y|y  =  0  =  0 

for  all  values  of  z  and  t,  and 

H(y,z=  ±k,t)  =  0  implying  that  h(y,z=  ±k,t)  =  0 


4 


These  conditions  hold  only  for  the  open  circuit  case  with  no  net  current  between 
conductors. 

Initially  the  total  velocity  and  induced  magnetic  intensity  are  zero,  which 
means  the  transient  velocity  and  the  transient  induced  magnetic  intensity  must  be 
the  negative  of  the  steady  solutions.  In  that  case 


u(y,z,t  =  0)  =  -Uss(y,z)  and 


h(y,z,t  =  0)  =  -Hss(y,z) 


MATHEMATICAL  SOLUTION 


To  automatically  satisfy  all  of  the  boundary  conditions  we  choose 
solutions  of  the  form 


u(y,z,t)  =  I  I  Umn(t)sin(rnTty)  cos(nnz/(2k)) 

m  =  0, 1 ,2..  n=  1,3.. 


h(y,z,t)  =  I  I  Hmn(t)cos(m7ty)  cos(nrrz/(2k)) 

m  =  0,l,2..  n=  1,3.. 


To  satisfy  the  partial  differential  equations  governing  the  solutions  for  u(y,z,t) 
and  H(y,z,t),  the  above  assumed  solutions  are  substituted  into  Eqs.  (7)  and  (8). 
There  results  a  pair  of  coupled  time  dependent  differential  equations  in  time 
given  by 


dUmn/dt  +  n2[m2  +  (n/(2k))2]Um  n  =  -MmrrH, 


PrndHmn/dt  +  Ti2[m2  +  (n/(2k))2]Hmn  =  MmnU, 


Equations  (11)  and  (12)  are  combined  to  yield  a  pair  of  second  order  time 
dependent  differential  equations  for  Umn  and  Hmn.  They  are  given  by 


Pmd2/dt2  +  (l+pm)[(mn)2  +  (nn/(2k))2]d/dt  + 


(Mmrt)2  +  [(mrt)2  +  (nn/(2k))2]2  „mn  =  0 


We  assume  the  solution  to  this  equation  is  Umn(t)  =  C  e_Pt  whence  P  must  be 
found  from  the  roots  to  the  equation 


PmP2  -  (l  +  Pm)aP  +  (Mmn)2  +  a2  =  0 


juuuuyiauk 


where  a  =  n2[m2  +  (n/(2k))2].  Because  the  two  solutions  of  Eq.  (14)  may  differ 
greatly  in  magnitude  we  choose  them  to  be  of  the  form8 


P,  =  [a(l  +Pm)  +  (1  -Pm)2a2  -  4Pm(Mmn)2]/(2Pm)  (15) 

and 

P2  =  [a  +  (Mmn)2]/(P[Pm)  (16) 

Thus  the  solutions  to  Eq.  (13)  are  given  in  exponential  form  by 

Umn(t)  =  C1e-p>t  +  C2e-P2t  (17) 

Hmn(t)  =  +  /3C2e~P2t  (18) 

where  a  =  (P(-a)/(Mmn)  and  ft  =  (P2-a)/(Mnn).  Utilizing  the 
aforementioned  initial  conditions  yields 

C|  +  Cj  =  -Ug„  (19) 

aC,  +  PC2  =  -H“„  (20) 

having  solutions  of 

C(  =  (U  £„-«„)/(*-/?)  (21) 

C2  =  (-aU«n  +  H«n)/(a-0)  (22) 


where  U“n  and  H”n  are  the  Fourier  coefficients  of  the  steady  solution1  when  it 
is  expressed  as  a  double  Fourier  series  of  the  forms  of  Eqs.  (9)  and  (10).  To 
reduce  the  amount  of  tedious  algebra  resulting  upon  the  invocation  of  the 
orthogonality  of  the  sinusoidal  functions  we  have  tried  to  express  the  steady 
state  results1  as  simply  as  possible.  These  solutions  are  given  by 


Uss(y,z)  =  1  Un(y)cos(n7tz/(2k))  (23) 

n  =  1,3,5... 

and 


Hss(y,z)  =  I  Hn(y)cos(nnz/(2k))  (24) 

n=  1,3,5... 


where  Un(y)  =  [g2(ei  ~  e2)/ci1  +  gi(e3_e4)/d2lF  » 

Hn(y)  =  f-gi(e1+e2)/d1  +  gi(e3  +  e4)/d2]F  , 


F  =  (4sin(nn/2))/((g[ +g2)n7t)  , 
gj  =  (M  +  N/M2  +  n¥/k2  )/2 
g2  =  (-M  +  VM2  +  n¥/k2)/2 
e,  =  e-g'(1“y) 
e2  =  e-g‘(1+y) 
e3  =  e-g2(,-y) 
e4  =  e~82(1+y) 
di  =  (1  -e)~2gl  and 
d2  =  (l-e)‘2g2. 

The  standard  use  of  orthogonality  yields  expressions  for  U“n  and  H^n  as  given 

by: 


1 

umn  =  -2  /  Un(y)  sin(mny)  dy  (25) 

0 

1 

Hmn  =  -  2  /  H„(y)  cos(mTty)  dy  (26) 

0 

The  indicated  integration  yields 

ss  _  (2cos(mn)mn(  - gf  +  gtg2  -  g2  -  m V2)) 

Umn  (g^g2  +  +  gjmV  +  m4;!4) 

and 

_  (2cos(mn)g1g2(g)  -g2)) 
mn  (g2g2  +  gfm2^2  +  g^m2^  +  mV*) 

The  substitutions  of  the  results  of  Eqs.  (27)  and  (28)  into  Eqs.  (21)  and  (22), 
and  the  quantities  defined  after  Eq.  (24)  and  considerable  amounts  of  algebraic 
manipulation  yields  C[  and  C2,  which  are  given  by 


Ci  =  (16cos(mn)sin((nn)/2)mPm(a2g^Op  -  a2g{g2Op  +  a2gjOp  +  a2m2Oprt2 

-  agfg2MOp  +  agfRa  -  2agf  +  ag^fMOp  -  agtg2Ra  +  2ag,g2 
+  agjRa  -  2agj  +  amVRa  -  2am 2n2  -  g^g2MRa  -  2gjM2m2rt2 

+  gtg|MRa  +  2gig2M2m2n2  -  2g^M2m2rr2  -  2M2m4Tr4))/(n(a2gjg20p 

2  2  2 

+  a2gim2OpTT2  +  a2g2m20prr2  +  a2m40p^  +  2ag|g20pRa 

-  4agjg2Pm  +  2ag2m2Opn2Ra  -  4ag2m2n2Pm  +  2ag2m2Oprt2Ra 

-  4ag2m2rT2Pm  +  2am40prr4Ra  -  4amVJPm  -  4g^M2m2n2Pm 
+  g2g^Ra  -  4g2M2m4rr4Pm  +  g2m2rt2Ra  -  4g2M2m4rt4Pm 

+  g^m2rt2Ra  -  4M2m6rt6Pm  +  m^Ra))  (27) 

and 

C2  =  (8cos(mn)sin((nn)/2)m(a2g^Op  -  2a2g^OpPm  -  a2g!g2Op 

2  2 

+  2ag!g2OpPm  +  a^Op  -  2a2g20pPm  +  a2m2Oprr2 

-  2a2m2Opn2Pm  +  2agjg2MOpPm  +  2ag^OpRa  -  2ag2PmRa 

-  2ag1g^MOpPm  -  2agjg2OpRa  +  2ag!g2PmRa  +  2ag2OpRa 

-  2agjPmRa  +  2am2Oprt2Ra  -  2am2rt2PmRa  +  2gfg2MPmRa  +  g,Ra 

-  2gtg^MPmRa  -  g^Ra  +  g^Ra  +  m2n2Ra))/(n(a2gfg§Op 
+  a2g^m2Opm2  +  a2g2m20p7t2  +  a2m40pn4  +  2ag^g|OpRa 

-  4agjg2Pm  +  2ag^m20pTT2Ra  -  4ag^m2n2Pm  +  2ag2m2Opn2Ra 

-  4ag2m2n2Pm  +  2am40pn4Ra  -  4am4n4Pm  -  4g^g|M2m2n2Pm 
+  gfg^Ra  -  4g^M2m47t4Pm  +  gjmVRa  -  4g^M2m4Ti4Pm 

+  gjm^Ra  -  4M2m6n6Pm  +  mV^Ra))  (28) 

For  documentation  we  give  the  expressions  for  crC)  and  /)C2  which  are  used  for 
the  calculation  of  the  induced  magnetic  intensity.  These  expressions  are 


(8cos(mrt)sin((nn)/2)(a3gjOp  -  2a3gfOpPm  -  a3gjg2Op 

1  2 

+  2a3g1g2OpPm  +  a3g^Op  -  2a3g|OpPm  +  a3m20pn2 

-  2a3m20pn2Pm  -  a2g2g2MOp  +  2a2g2g2MOpPm  +  2a2gfOpRa 

-  2a2g2Op  -  2a2g2PmRa  +  4a2gfPm  +  a2gjg^MOp 

-  2a2g1g2MOpPm  -  2a2g1g20pRa  +  2a2g!g20p  +  2a2g!g2PmRa 

-  4a2g]g2Pm  +  2a2g2OpRa  -  2a2g^Op  -  2a2g^PmRa  +  4a2g2Pm 
+  2a2m2Opn2Ra  -  2a2m2Opn2  -  2a2m27t2PmRa  +  4a2m2it2Pm 

-  2ag|g2MOpRa  +  2agjg2MPmRa  -  2ag2M2m2Opn2 

+  4ag2M2m2n2Pm  +  ag^Ra  -  2agfRa  +  2agtg^MOpRa 

-  2agtg^MPmRa  +  2ag1g2M2m20pn2  -4ag1g2M2m2rt2Pm  -  ag,g2Ra 
+  2agjg2Ra  -  2ag|M2m2Opn2  +  4ag|M2m2n2Pm  +  ag^Ra  -  2ag|Ra 

-  2aM2m4Opn4  +  4aM2m4n4Pm  +  am¥Ra  -  2am2n2Ra 

-  g2g2MRa  -  2g2M2m2rt2Ra  +  gjg^MRa  +  2g1g2M2m2n2Ra 

-  2g^M2m2n2Ra  -  2M2m4tr4Ra))/(Mnn(a2gfg|Op  +  a2gfm2OpTi2 
+  a2g^m2Oprr2  +  a2m40pn4  +  2agfg|0pRa  -  4ag2g^Pm 

+  2ag2m2OpTi2Ra  —  4ag2m2rt2Pm  +  2ag2m2Opn2Ra  —  4ag2m2n2Pm 
+  2am4Optt4Ra  -  4am4TT4Pm  -  4g2g2M2m2n2Pm  +  gfg2Ra 

-  4g2M2m47J4Pm  +  g^mVRa  -  4g^M2m4n4Pm  +  g^m27i2Ra 

-  4M2m6Ti6Pm  +  mV^Ra))  0 

=  (8cos(mn)sin((nrr)/2)(  -  a3g  jOp  +  2a3g^OpPm  +  a3g(g20p 

-  2a3g1g2OpPm  -  a3g20p  +  2a3g^OpPm  -  a3m20prr2 

+  2a3m2Opn2Pm  -  2a2g?g2MOpPm  -  2a2gfOpRa  +  2a2gfOp 
+  2a2gfPmRa  -  4a2g2Pm  +  2a2gjg2MOpPm  +  2a2gjg2OpRa 


-  2a2g1g2Op  -  2a2g1g2PmRa  +  4a2g,g2Pm  -  2a2g2<I)pRa  +  2a2g20p 
+  2a2g2PmRa  -  4a2g|Pm  -  2a2m2Opn2Ra  +  2a2m2Opn2 

+  2a2m2tt2PmRa  -  4a2m2n2Pm  -  2ag2g2MPmRa  +  4ag2g2MPm 
+  2ag^M2m2Opn2  -  4ag^M2m2rt2Pm  -  ag2Ra  +  2agjRa 
+  2agig2MPmRa  -  4ag[g2MPm  -  2ag1g2M2m20pTt2 
+  4ag1g2M2m2rr2Pm  +  agig2Ra  -  2ag[g2Ra  +  2ag2M2m2Oprt2 

-  4ag2M2m2n2Pm  -  ag2Ra  +  2ag^Ra  +  2aM2m4Oprt4 

-  4aM2m4Tr4Pm  -  am2n2Ra  +  2am2n2Ra  +  4g2g2M3m2rr2Pm 

+  2g2M2m2n2Ra  -  4g1g|M3m2rt2Pm  -  2g1g2M2m2rt2Ra  +  2g2M2rt2Ra 
+  2M2m4Tr4Ra))/(Mnrt(a2g2g2Op  +  a2gjm2Opn2  +  a2g2m2Optt2 
+  a 2m4Opn4  +  2agjg20pRa  -  4ag2g2Pm  +  2ag2m2Opn2Ra 

-  4ag2m2n2Pm  +  2ag2m2Oprr2Ra  -  4ag2m2n2Pm  +  2am4Opn4Ra 

-  4am4n4Pm  -  4g^M2m2n2Pm  +  g2g^Ra  -  4g^M2m4n4Pm 
+  g^mVRa  -  4g^M2m4n4Pm  +  g^mVRa  -  4M2m6n6Pm 

+  mVRa))  (30) 

The  coefficients  of  the  time  functions  in  the  exponentials  are  given  by 

PI  =  (m'/64M2k4m2Pm  -  16k4m4n2Pm+  32k4m4n2Pm  -  16k4m4n2 
-  8k2m2n2n2Pm  +  16k2m2n2tt2Pm  -  8k2m2n2n2  -  n47i2Pm 


+  2n4n2Pm  -  n4n2)  +  4k2m2nPm  +  4k2m2n  +  n^Pm 
+  n  ^/(SkPm) 


(?D 


P2  =  (2n(4M2k2m2  +  4k2m2  +  n2))/(i'/(64M2k4m2Pm  -  16k4m4n2Pm 


+  32k4m4n2Pm  -  16k4mn2  -  8k2m2n2n2Pm  +  16k2m2n2n2Pm 


-  8k2m2n2n2  -  n4n2Pm+  2n4n2Pm  -  n4n2)  +  4k2m2rcPm 
+  4k2m2n  +  n^Pm  +  n2^)  (32) 


M  =  Hartmann  number,  Op  =  1  +  Pm,  Om  =  1  -  Pm  and  Ra  = 
'J Oma2  -  4Pm(M7im)2 


Thus  the  transient  parts  of  the  velocity  and  magnetic  intensity  mathematical 
expressions  are  determined  by  Eqs.  (9)  and  (10)  and  the  auxiliary  Eqs.  (17),  (18), 
and  (27)  through  (32).  In  addition  the  above  given  definitions  of  gj ,  g2.  Op, 

Om,  and  Ra  must  be  used. 

LIMITING  CASE  SOLUTIONS 

When  the  Hartmann  number  approaches  zero  the  solution  given  above 
approaches  the  nonmagnetohydrodynamic  limit  of  zero  induced  magnetic 
intensity  and  a  purely  hydrodynamic  transient  fluid  dynamic  velocity 
distribution.  The  complete  nonmagnetohydrodynamic  solution  (including  both 
the  steady  and  the  transient  parts)  is  given  by 

OO 

U(y,z,t)  =  4  I  [sin(mt/2)cos(nTtz/(2k))/[7m(l  -  e ~ nTt/k)]] 
n=  1,3.. 

[enn(y-  l)/(2k)  _  e-nrt(y+  l)/(2k)] 


00  m  cos(mn)sin(n/2)cos(nnz/(2k))e  ~  tff2lm2  +  (n/(2k))2] 
m=  1,2,3..  n2n  [m2  +  (n/(2k))2  (33) 

In  most  laboratory  experiments  Pm  is  of  the  order  of  2x  10~6  and  the 
expressions  for  Pt  and  P2  simplify  to 

Pj  =  n2[m2  +  n2/(4k2)]2/Pm  (34) 


and 


P2  =  1  +  M2m2/[m2  +  n2/(4k2)] 


(35) 


where  it  is  clear  that  the  P2  terms  contained  in  the  second  time  function  of  Eqs. 
(9)  and  (10)  do  not  aid  in  the  convergence  of  the  series,  which  is  made  up  of 


double  summations.  On  the  other  hand  the  Pj  terms,  in  addition  to  being  very 
large,  increase  in  size  with  the  square  of  the  indices  m  and  n.  From  the 
inspection  of  these  variables  it  is  expected  that  the  first  term  of  the  exponential 
solutions  not  only  converges  faster  with  index  but  also  vanishes  faster  than  the 
second  term. 

Exponential  diminishing  time  solutions  were  assumed  for  the  temporal 
variations  as  a  matter  of  convenience.  However,  inspection  of  the  radical 
contained  in  the  solution  of  Pj  indicates  the  possibility  of  complex  roots.  It  is 
clear  that  oscillatory  solutions  exist  when 

M  >  (1  -  Pm)[m2  +  n2/(4k2)]/(2mrtPm)  (34) 

requiring  extremely  large  Hartmann  numbers  in  the  terrestrial  laboratory, 
assuming  high  temperature  liquid  metal  superconductors  are  not  available.  If  the 
magnetic  Prandtl  number  can  be  increased,  then  some  of  the  early  terms  in  the 
summations  may  be  damped  sinusoids. 

RESULTS  AND  DISCUSSION 

We  present  a  few  calculations  of  velocity  and  induced  magnetic  intensity 
transients,  roughly  covering  the  area  of  the  channel.  The  calculations  are  done 
in  double  precision  and  require  the  summation  of  50,000  terms  for  each  point  in 
space9  (m  and  n  ranging  up  to  1000  and  50  respectively).  By  the  use  of 
minimum-maximum  method10  only  400  terms  are  necessary  for  engineering 
accuracy.  The  optimization  of  the  method  for  double  series  summations  might 
further  reduce  the  needed  number  of  terms.  To  focus  our  attention  upon  the 
temporal  variation  of  the  velocity  and  the  induced  magnetic  intensity,  both 
quantities  are  normalized  by  dividing  them  by  their  final  or  steady  state  value. 
The  steady  state  values  are  the  same  as  those  presented  in  an  earlier  paper.1  By 
this  technique  the  transients  at  different  points  in  the  channels  can  be  compared 
most  easily. 

TRANSIENTS 

In  Figs.  2,  4,  6  and  8  are  shown  normalized  dimensionless  velocity 
transients  for  an  aspect  ratio  of  unity  a  Hartmann  number  of  5  and  a  magnetic 
Prandtl  number  of  2.54  x  10 -6.  In  each  case  the  velocity  was  normalized  with 
respect  to  its  steady  state  value.  In  each  figure  there  are  curves  of  normalized 
dimensionless  velocity  as  a  function  of  dimensionless  time  at  y  =  0.25,  0.5  and 
0.75.  A  different  value  of  z  is  used  in  each  of  the  figures.  In  the  case  where  z 
=  0  (see  Fig.  2),  the  transients  are  almost  identical  at  all  values  of  y  presented. 
Inspection  of  Figs.  2,  4,  6  and  8  indicates  the  differences  as  the  transients 
increase  slightly  when  z  increases  from  the  channel  center  to  the  insulating  walls 
and  the  risetime  of  the  transient  also  increases  slightly.  The  risetime  for  the 
velocity  in  the  cases  presented  here  is  roughly  0. 1  units  of  dimensionless  time. 
Also,  for  the  velocity  transients  the  exponential  solution  having  the  small  time 
constant  has  small  amplitude  compared  to  the  large  time  constant  exponential 
solution.  That  is,  C(  is  less  than  C2. 


These  curves  are  presented  in  Figs.  3,  5,  7  and  9  for  the  same  parameters 
as  used  in  the  corresponding  velocity  figures.  In  each  case  the  induced  magnetic 
field  was  normalized  with  respect  to  its  steady  state  value.  In  all  cases  the 
normalized  dimensionless  induced  magnetic  intensity  rises  rapidly  from  zero  to 
an  apparent  initial  value  (see  Fig.  3)  and  subsequently  undergoes  less  rapid 
changes.  The  expansion  of  the  dimensionless  time  scale  shows  that  the 
normalized  dimensionless  induced  magnetic  intensity  actually  begins  at  a  value 
of  zero  (as  it  should)  and  increases  very  rapidly  at  first.  This  is  the  small  time 
constant  exponential  part  of  the  solution  and  as  it  fades  the  long  time  constant 
exponential  part  of  the  solution  takes  over.  The  magnetic  intensity  initially 
builds  up  faster  at  stations  nearer  the  moving  plate.  We  feel  that  the  fast  part  of 
these  transients  is  due  to  the  passing  of  Alven  waves  while  the  slow  part  is  due 
to  viscous  and  magnetic  diffusion.  The  curves  are  similar  for  the  different  values 
of  z  and  show  a  slight  tendency  to  have  a  slower  approach  to  the  steady  state  as 
z  increases. 


CONCLUSION 

A  formal  solution  to  the  transient  magnetohydrodynamic  Couette  problem 
has  been  obtained  and  we  have  presented  but  few  of  the  situations  that  could  be 
examined.  We  have  disclosed  the  presence  of  fast  and  slow  transients  and  the 
possibility  of  oscillatory  behavior  in  both  the  fluid  velocity  and  the  induced 
magnetic  field.  In  future  work  we  plan  to  investigate  these  phenomena  and  start¬ 
up  losses. 
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Fig.  3.  Magnetic  intensity  transients  at  z  =  0,  y  =  0.25,  0.5  and  for  k  *  1  and 
a  Hartmann  number  of  5. 

0.75 


Fig.  4.  Velocity  transients  at  z  ■  0.25,  y  =  0.25,  0.5  and  0.75  for  k  *  1  and 
a  Hartmann  number  of  5. 
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DIMENSIONLESS  TIME,  t 

Fig.  5.  Magnetic  intensity  transients  at  z  =  0.25,  y  =  0.25,  0.5  and  0.75, 
for  k  *  1  and  a  Hartmann  number  of  5. 


DIMENSIONLESS  TIME,  t 

Fig.  6.  Velocity  transients  at  z  *  0.5,  y  =  0.25,  0.5  and  0.75  for  k  =•  1  and 
a  Hartmann  number  of  5. 
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